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Express Label EL6271 3 1 377US ! 



Date of Deposit: June 20, 2000 

DETERMINING THE EQUIVALENCE OF TWO SETS OF SIMULTANEOUS 

LINEAR ALGEBRAIC EQUATIONS 



The present invention relates to a computer implementable method, and in 
particular, to a method and apparatus for determining whether two sets of simultaneous 
linear algebraic equations are equivalent. 



In many applications, the need arises to solve one or more systems of 
simultaneous linear algebraic equations (SLAEs) whose coefficient matrices comprise 
only numerical elements. Such applications include engineering and simulation computer 
codes. Solutions of the SLAE are typically obtained by using the well-known Gaussian 
elimination method. Therefore, prior methods typically would solve two such SLAE 
systems S\ and S 2 , and compare their solutions. However, such methods may not always 
work if one or both of the SLAEs are ill-conditioned and/or the numerical precision used 
in computations is not high enough. 

Furthermore, such methods are generally not adapted to solving a set of SLAEs 
whose coefficient matrix elements are algebraic expressions, and for which the solution 
will, in general, be in algebraic form. 



It is an object of the present invention to provide a method of determining whether 
two sets of simultaneous linear algebraic equations are equivalent. 

The invention provides a computer implemented method for determining the 
equivalence of two sets of simultaneous linear algebraic equations (SLAEs), each of the 
sets comprising two or more algebraic equations. The method comprising the steps of: 

reducing each SLAE to a standard form; and 
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comparing the SLAEs to determine whether equivalence exists. 

The invention further provides a computer implemented method of determining the 
equivalence of a first and a second set of simultaneous linear algebraic equations 
(SLAEs), with the method comprising the steps of: 

iteratively eliminating unknowns from each of the sets of SLAEs to place each 
SLAE in a two-part standard form; and 

forming a product of a part of one standard form equation with a part of another part 
of another standard form equation; 

forming a product of the other part of standard form equation with the other part of 
another standard form equation; and 

comparing the respective products for mathematical equivalence. 

There is further provided a computer implemented method of determining the 
equivalence of a first and a second set of simultaneous linear algebraic equations, each of 
the equations being of a form: 

e±ixi + e i2 ^2 + e i3 x 3 + . . . + e in x n = b ± 
wherein Xj are unknowns, e±j are coefficients, and b± are quantities. The coefficients 
and quantities are known algebraic expressions. The method comprising the steps of: 

iteratively eliminating the unknowns from each of the sets of simultaneous linear 
algebraic equations until each of the equations are in the form: 



wherein la and r± are algebraic expressions, and £={1;2} indicate one of the sets that the 
equation is derived from; and 

comparing, for each of the unknowns, the products (l±i)i*(ri) 2 and 
da) 2* i, wherein the first and the second set of simultaneous linear algebraic 
equations are equivalent if the products match for all the unknowns. 

The invention further discloses a computational apparatus for determining the 
equivalence of a first and a second set of simultaneous linear algebraic equations, each of 
the equations being in the form: 

eiiXi + e i2 x 2 + e i3 x 3 + . . . + e in x n = b± 



da)k Xi = (r ± ) k 
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wherein Xj are unknowns, e±j are coefficients, and hi are quantities, the coefficients and 
quantities being known algebraic expressions. The apparatus comprising: 

means for iteratively eliminating the unknowns from each of the sets of 
simultaneous linear algebraic equations until each of the equations are in the form: 

(la) k x ± = (r±) k 

wherein la and r± are algebraic expressions, and £={1;2} indicate one of the sets that the 
equation is derived from; and 

means for comparing, for each of the unknowns, the products ( 1±±) i* (r±) 2 and 
da) 2* (r±) i, wherein the first and the second set of simultaneous linear algebraic 
equations are equivalent if the products match for all the unknowns. 

The invention yet further discloses a computer program product carried by a storage 
medium for determining the equivalence of a first and a second set of simultaneous linear 
algebraic equations, each of the equations being of a form: 

e± 2 x x + e i2 x 2 + e i3 x 2 + . . . + e in x n - b± 
wherein Xj are unknowns, G±j are coefficients, and b± are quantities, the coefficients and 
quantities being known algebraic expressions. The computer program product 
comprising: 

a program element for iteratively eliminating the unknowns from each of the sets of 
simultaneous linear algebraic equations until each of the equations are in the form: 

da) k *i = (r±) k 

wherein la and r± are algebraic expressions, and k={\\2) indicate one of the sets that the 
equation is derived from; and 

a program element for comparing, for each of the unknowns, the products 
da) 1* 2 and ( la) 2 * (*i) i, wherein the first and the second set of simultaneous 
linear algebraic equations are equivalent if the products match for all the unknowns. 

Preferably, the method further includes recasting the algebraic expressions into a 
form of one or more token pairs arranged sequentially in a string, each of the token pair 
comprising an operator followed by an operand; and reducing the strings in accordance 
with a set of predetermined simplifying rules to obtain reduced expressions. Eliminating 
the unknowns from each of the sets of simultaneous linear algebraic equations is 
performed on the reduced strings in accordance with a set of predetermined operations. 
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Furthermore, the simplifying rules can comprise the steps of arranging token pairs 
into subgroups, arranging operand tokens in an arranged subgroup in order, reducing the 
ordered operands by consolidating one or more constants and eliminating variables of 
opposite effect to form reduced subgroups, and consolidating one or more multiple 
instances of similar subgroups, to produce a reduced string. 



A preferred embodiment of the present invention will now be described with 
reference to the drawings in which: 

Fig. 1 is a schematic block diagram of a conventional general-purpose computer 
system upon which the embodiment of the invention may be practised; and 

Fig. 2 is a flow diagram of a method of determining whether two sets of 
simultaneous linear algebraic equations are equivalent. 



A general-purpose computer system 100, upon which the preferred embodiment 
of the invention may be practised, is shown in Fig. 1. The computer system 100 will first 
be described, followed more particularly by a description of a method of determining 
whether two sets of simultaneous linear algebraic equations are equivalent. 

This method may be implemented as software, such as an application program 
executing within the computer system 100. In particular, the steps of the method of 
determining whether two sets of simultaneous linear algebraic equations are equivalent, 
are effected by instructions in the software that are carried out by the computer system 
100. The software may be stored in a computer readable medium, including the storage 
devices described below, for example. The software is loaded into the computer system 
100 from the computer readable medium, and then executed by the computer system 100. 
A computer readable medium having such software or computer program recorded on it 
is a computer program product. The use of the computer program product in the 
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computer preferably effects an advantageous apparatus for determining whether two sets 
of simultaneous linear algebraic equations are equivalent, in accordance with the 
embodiments of the invention. 

The computer system 100 comprises a computer module 101, input devices such 
as a keyboard 102 and mouse 103, and output devices including a printer 115 and a 
display device 114. The computer module 101 typically includes at least one processor 
unit 105, a memory unit 106, for example formed from semiconductor random access 
memory (RAM) and read only memory (ROM), input/output (I/O) interfaces including a 
video interface 107, an I/O interface for the printer device 115 and an I/O interface 113 
for the keyboard 102 and mouse 103. A storage device 109 is provided and typically 
includes a hard disk drive 110 and a floppy disk drive 111. A CD-ROM drive (not 
illustrated) may be provided as a non- volatile source of data. The components 105 to 113 
of the computer module 101, typically communicate via an interconnected bus 104 and in 
a manner which results in a conventional mode of operation of the computer system 100 
known to those in the relevant art. 

Typically, the application program of the preferred embodiment is resident on the 
hard disk drive 110, and read and controlled in its execution by the processor 105. 
Intermediate storage of the program may be accomplished using the semiconductor 
memory 106, possibly in concert with the hard disk drive 110. hi some instances, the 
application program may be supplied to the user encoded on a CD-ROM or floppy disk 
and read via a CD-ROM drive (not illustrated) or floppy disk drive 111, or alternatively 
may be read by the user from the network (not illustrated) via the modem device (not 
illustrated). Still further, the software can also be loaded into the computer system 100 
from other computer readable medium including magnetic tape, a ROM or integrated 
circuit, a magneto -optical disk, a radio or infra-red transmission channel between the 
computer module 101 and another device, a computer readable card such as a PCMCIA 
card, and the Internet and Intranets including e-mail transmissions and information 
recorded on websites and the like. The foregoing is merely exemplary of relevant 
computer readable mediums. Other computer readable mediums may be practiced 
without departing from the scope and spirit of the invention. 
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Having described the hardware environment of the invention, the method of 
determining whether two sets of simultaneous linear algebraic equations are equivalent 
will now be described. 

Broad Outline of Method 

Let S represent a system of simultaneous linear algebraic equations (SLAEs) as is 
given by the following: 

e x x x\ + e x2 x 2 + e\yxz + ...+ e\„x„ = b\ 
e 2 \x\ + e 22 x 2 + e 23 x 3 + -..+ e 2n x n = b 2 

e n ix\ + e n2 x 2 + e„3^3 + ■...+ = 6„ 

where w-unknowns {xx&jci, are related by n equations, and coefficients e,y (with 

z=l,2,...,« and 7= 1,2,. are known algebraic expressions, as are the right-hand side 
quantities Z>„ i = 1,2,...,«. 

The method of determining whether two such systems S\ and S2 are equivalent - 
that is, their respective solutions are identical to each other - broadly has two parts, 
namely: 

(1) The reduction of each system of SLAEs S into a standard form of the type 

l\\X\ — r\ 
hi*i — r 2 
hi*! = n 

where /„- and r,* are algebraic expressions; and 

(2) Comparison of two sets of SLAEs in their standard form. 

It is assumed that the coefficients ey and the quantities 6, of the SLAEs S\ and S 2 
have no division operators. Undesirable division operators can be eliminated from the 
SLAEs S\ and S 2 by multiplying the affected equations by appropriate factors. This is 
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done to reduce the complexity of handling operands associated with the division operator, 
which is not a commutative operator. 

Reduced Expression 

The coefficients e^ and the quantities 6, may be written as expressions, wherein the 
terms in the coefficients ey and the quantities bi may include constants and variables. In 
the preferred embodiment, to facilitate comparisons between two expressions, the concept 
of a reduced form of an expression, as described below, has been used. The reduced 
expression is the canonical form to which expressions are converted. 

It is apriori assumed that the expression to be converted is syntactically correct and 
does not contain any blanks. In the preferred embodiment, variables are limited in their 
construction to lower-case alphabets, underscore character, and digits, except that a 
variable may not start with a digit or end with an underscore. If these construction rules 
are not met, then the affected variables may be mapped (aliased) to alternative, but 
distinct, variables obeying the construction rules, and these new variables used instead. 

A convention adopted for the present embodiment is that variables in the 
coefficients e^ and the quantities b t raised to a positive integer power are written out as 
multiplications of the variables. Thus, for example: 

a 11 becomes a* a* . . . *a, where a appears n times in the product. 

To convert a given expression into a reduced expression, the expression firstly is 
put in the following form: 

<unitary operator><operand><operator><operand> <operatorxoperand> 

where the unitary operator is either + (plus) or - (minus), and each operator is one of + 
(plus), - (minus), or * (multiplication). In the event that an expression does not 
commence with a unitary operator, a unitary operator + (plus) is inserted at the start of the 
expression. For example: 

a+b*c-d becomes +a+b*c-d 

Note, in particular, the absence of brackets. Brackets, if present, in the 
expression, must be removed by carrying out the necessary operations needed to remove 
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them, such as multiplying two parenthesized factors, discarding superfluous brackets, etc. 
to bring a given expression into the above form. 

Next, all + (plus) operators are substituted with the string +1* so that + becomes 
+ 1*. Similarly, all - (minus) operators are substituted with the string -1* so that - 
becomes -1*. Thus, for example: 

+a becomes + 1 * a 



Finally, the operands, which are constants (including the 'l's introduced in the 
previous step) are converted into an e- format as follows: 

".[unsigned number] e[e-sign] [unsigned exponent] 1 ' 

where: [unsigned number] is a n-digited number comprising only digits and n is a 
prefixed integer greater than 0; 

[e-sign] is the sign of the exponent and is one of > for plus or < for minus; and 
[unsigned exponent] is a m-digited number comprising only digits and m is a 
prefixed integer greater than 0. 

Thus, for example: 

25 = 0.25*10 2 becomes . 250000e>02 

and 

0.025 = 0.25*10 -1 becomes .250000e<01 

where it is assumed n=6 and m=2. It is noted that any constant will be represented by a 
string of constant length m+n+3 characters in the e- format. Here e[e-sign] [unsigned 
exponent] represents the quantity 10 raised to the power [e-sign] [unsigned exponent], 
which must be multiplied to the number represented by .[unsigned number] to get the 
actual constant. 



and 



-a*b 



becomes 



-l*a*b 
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Now, the expression will contain at least one operand which is a constant. Each 
expression will have one or more terms, where each term has the following form: 

<unitary operator><operand><*xoperand> <*><operand> 

where the unitary operator is either + (plus) or - (minus), and between two consecutive 
operands is the multiplication operator *. After the terms are identified, the [e-sign] of 
each constant is restored from < or > to - or + respectively. 

In each term the operands are sorted (rearranged) in ascending order according to 
their ASCII (American Standard Code for Information Interchange) value. This does not 
affect the term since the multiplication operator is a commutative operator, so the 
exchange of operands is completely permissible. The operands, which are constants, will 
all bunch up at the beginning of the terms where they can be easily identified and replaced 
by a single constant. Thus, for example: 

+. 10 000 0e+01*a*b* . 500000e+00 
after arranging the operands in ascending order becomes 

+ . 10 0 00 0e+01* . 5 0 000 0e+0 0*a*b 
and after consolidating the constants the term becomes 

+ . 500000e+00*a*b 

At this stage a term will have the following form: 
<unitary operator><constant><*><operand> <*><operand> 

where each operand is a variable, whose ASCII value is not lower than that of its 
preceding operand, if any. This is the reduced form of a term. In the reduced form, the 
non-constant part of a term is called a variable-group. For example, if the term in the 
reduced form is " + . 250000e+01*a*a*b n , then its variable-group is " *a*a*b". 

In an expression, all those terms whose variable-groups match, are combined by 
modifying the constant in one of the terms, and eliminating all other terms with identical 
variable-group. 
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Finally, the reduced terms in the expression are rearranged in an ascending order 
according to the ASCII value of their respective variable-group. In this final form, the 
expression is said to be in its reduced form. Note, in particular, that no two terms in a 
reduced expression will have the same variable-group. 

Method of Determining Equivalence 

Referring to Fig. 2, a method 200 of determining whether two such systems S\ and 
S 2 are equivalent is shown. Starting in step 240, all the coefficients ey and the quantities 
bt are converted into their respective reduced form (as discussed above). 

In steps 250 to 280, the Gaussian elimination and back substitution method 
(adapted to avoid divisions) is used to bring the SLAEs S\ and S 2 into a standard form. 

In step 250 a counter k is set to 1. Step 252 follows, where the variable x k is 
eliminated from the y-th equations, j = (k+l),...,w, to get a &'th derived system. In 
particular, with counter k equal to 1, the variable xi is eliminated from the y-th equations, 
j = 2,3,. to get a first derived system defined as: 

e\\x\+ e\ 2 x 2 + ei3*3+...+ e Xn x n = b\ 
l e 2 2*2 + ^23*3 + ... + l e 2n x„ = ] b 2 

x e n2 x 2 + x e n3 x 3 + - + X e nn x n = x b n 

where the new coefficients x e jk of the first derived system are given by: 

X e jk = ejk e x \ - e\ k ej\ ; and 

x bj = bjen - b\ e^ for (j,k) = 2,...,n. 

In a case where the coefficient e\\ =0, then the first equation of the system S is 
interchanged with any other equation m of the system S for which its coefficient e ]m is 
non-zero. If no such equation m can be found, then the SLAEs are singular, and the 
method 200, and in particular step 252, is interrupted by following the line 262 to step 
270, where the method 200 is terminated with an appropriate error message. 



1 
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In step 260 it is determined whether the counter k is equal to n-l, where n is the 
number of unknowns. If this is not so, a sub-system is defined in step 255 from the Ar'th 
derived system. For example, with counter k equal to 1, the sub-system derived from the 
first derived system is as follows: 

x e 22 x 2 + X e 23 X3 + ... + X e 2n *n = l b 2 
X e n2 x 2 + x e n3 *3 + + X e m *n = l b n 

This sub-system is a set of (n - 1) SLAEs, having (n - 1) unknowns {x^i,...^} . 
After incrementing the counter k in step 258, the steps of reduction 252 to 260 are now 
repeated on the sub-systems, until the system S is reduced to a (n -l)-th derived system as 
follows: 

e u x x + e X2 x 2 + 6fi 3 x 3 +... + e in x n = b\ 
x e 22 x 2 + x e 23 x 3 + . . . + 1 e 2n x n = X b 2 



nA e x = nA h 



wherein the diagonal coefficients y— are all nonzero, and where 

/ _ /-i i-\ /-i i-i 

ejk - ejk en - eut eju 

% = lA bj lA e u - M 6/ lA e jh for / = - 1; (j,k) = I +!,...,«, 



and 



ejk = ejk. 



This completes the Gaussian elimination phase of the process. Note the absence 
of any division in the entire process. The counter k is now equal to w-1 and the method 
therefore continues to step 280 where back substitution is performed, again without any 
division. Therefore, instead of calculating the unknown jc,-, the product is calculated, 
where each of the n unknowns x t is expressed in the form of a ratio x/ = r i /l ii with r, a 
numerator and /,,* a denominator. With i = n, we have, 

^nn-^n u n 

so that 

Un = n ' X e nn and r n = n ~ x b n . 
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For i = n - 1, the (n - l)-th equation is multilpied by the denominator l nn to obtain 

or 

/ AI_2 / o _ «-2l / n-2 _ 

so that 

= inn n 2 e n -\ yTl -\ and r„.i = " 2 b n -\ l nn - nJ1 e n -\si r n 
For / = n - 2, we multiply the {n - 2)-th equation by the denominator l n .\,n-\ and 

obtain 

n 3 ^/i-2 -*/i-2 + ln-\,n-\ " 3 ^n-2 t n-\ Xn-\ + ln-\ji-\ " 3 ^n-2 i n *n = " 3 &/i-2 

or 

e n -2,n-2 Xn-2 ~ On-2 ln-\,n-\ " ^n-\ y nA ^n-2,n ?n ~ ^n-2 y n-\ ^n-l 

so that 

/«-2,/i-2 = //i-u/i-i /i_3 ^«-2,/ I -2 and r„_ 2 = n ' 3 b n . 2 In-i^-i - rt " 2 ^-i,/i-i n ' 3 e n , 2ytl r n - n ~ 3 e n - 2s n-\ r n .\ 

It can be shown that for any / = 1 5 2,...,« - 1, the result will be 

la = l ' X eu and n = lA bi - R in r n - R in .i r nA - ... - R u + X r i+ \ 

with 

lnn = n ' X e nn and r n = n ~ l b n 

where 

Rij = lA eij for/ = «,...,(/ +1) and / = 1,2,.. .,n - 1. 

Note that since Ijj is a factor of Ry will be free of any divisions. However, it 

is noted that there is no step in the back substitution step 280 where factors common to /„ 
and ri have been eliminated. 

After completing steps 240 to 280 for each of the two SLAEs systems S\ and S 2 , 
string arrays (/,-,•) i and for system S\ and (ln) 2 and (n) 2 for system S 2 have been 
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produced. In principle, to show that the solutions of the two systems S\ and S2 are 
equivalent, it would suffice if their respective string arrays /,, and n where shown to 
match. However, this can not always be done, since it is generally not possible to 
eliminate their common factors completely by presently known methods. It must 
therefore be assumed that there may be uneliminated common factors present. However, 
it is clear that mathematically 



in which form a comparison may be performed. Therefore, step 290 calculates 
expressions (/,-,-) 1 * (V/) 2 and (//;) 2 * for each i = !,..,«. If all the expressions (/,-,-) 1 * (r,) 2 
and (7,7)2 * (n)\ have been consistently reduced to their reduced form, then a step 300 
performs a simple string comparison of (/,-,•) 1 * (r,) 2 with (7 i7 ) 2 * (r,-)i. A decision step 310 
determines whether matches were found for all i = 1,..,«. If the answer is Yes, then 
equivalence of systems Si and S 2 is reported in step 312. Alternatively, non-equivalence 
is reported in step 315. 



An example of performing the method 200 to determine whether two systems S\ 
and #2 are equivalent, will now be described. C and C++ programing language notations 
will be used. In this notation, the coefficients e^ and the quantities bi are denoted as 
e[i-l][j-l] andb[i-l] respectively. 

To understand the example given below, reference to the following pseudo-code 
fragment will be helpful. The variables e [ ] [ ] and b [ ] are assumed to have the 
datatype algebraic "expression", which inter alia will implement the operators '+' (plus), 
'-' (minus), and (multiplication) operators on such expressions. The class 
"expression" also has a method which can convert an algebraic expression into its 
reduced form. 

// Gaussian elimination 

// e[] [] and b[] are of type Expression. 

for (i = 0; i < n-1; i + + ) { // Index for the derived system. 



or equivalently, 



Example 
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// Comment 1 

// If e[i] [i] = 0, exchange this row with another below it (say 

// the k-th row, k > i) such that e[i] [k] 1= 0. If no such k is 

5 // found, exit with the message that the matrix e is singular. 

// The code to do this is not shown here. 

for (j = i+l; j < n; j++) { 
for (k = i+l; k < n; k++) { 

10 // Multiply i-th row with e[j] [i] . 

// Multiply j-th row with e[i] [i] . 
// Subtract i-th row from j-th row. 

e[j][k] = e[j] [k]*e[i] [i] - e [i] [k] *e [ j ] [i] ; 

} 

15 b[j] = b[j] *e[i] [i] - b[i] *e[j] [i] ; 

} 

// Zero lower triangle coefficients 
for (k = 0; k < i ; k++) e[i] [k] = "0"; 

20 } 

/ / Back-substitute . 
i = n; 

// At the end of the following while loop, e[i-l] [i-1] will 
25 // contain l ±i 

// and b[i-l] will contain r±. The solution will be x± = la/ r ± „ 

while ( i - - ) { 
j = n; 

30 while (j--) btj] - b[j] *e[i] [i] - b[i] *e[j] [i] ; 

for(k =0; k < n; k++) { 
for (j = k; j < n ; j++) { 
e[k] [j] *= e[i] [i] ; 

} 

35 } 
} 



Now, let system S\ be the set of equations: 
40 ax\ + X2 + X3 = a + 2 

X\ + X2 + X3 = 3 
•^1 ^2 ~ X3 — 1 

and let system ^2 be the set of equations 

ax\ + 2x2 = a + 2 
45 2x\ + 2x 2 = 4 

X2 — X3 =0 
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That is, each set consists of three equations. 



Considering system S\ first, the coefficients e,y and the quantities b t may be 
written as follows: 

e[0][0] = a 
e[0][l] = 1 
e[0][2] = 1 
b[0] = a+2 
e[l][0] = l 
e[l][l] = l 
e[l][2] = 1 
b[l] = 3 
e[2][0] = 1 
e[2][l] = 1 
e[2][2] = -l 
b[2] = 1 

Performing step 240 in system Sj, all the terms in the coefficients e,y and the 
quantities 6, are converted by the computer program performing method 200 into the 
reduced form, with the text variable to which a pseudocode variable refers to at different 
stages of computation noted on the right hand side, as follows: 



Reduced Form Variables 

e[0][0] = +.10000e+01*a °e n = e u 

e[0][l] = +.10000e+01 °e, 2 = ei2 

e[0][2] = +.10000e+01 °e 13 = e^^ 
b[0] = +.10000e+01*a+.20000e+01 = b x 

e[l][0] = +.10000e+01 °e 2 \ = e 2 \ 

e[l][l] = +.10000e+01 °e 22 = e 22 

e[ 1 ] [2] = +. 1 0000e+0 1 °e 23 = e 23 

b[l] = +.30000e+01 °b 2 = 6 2 

e[2][0] = +.10000e+01 °e 3 , = e 3 . 
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e[2] [ 1 ] = +. 1 0000e+0 1 °e 32 = e 32 

e[2][2] = -.10000e+01 °e 33 = e 33 

b[2] = +.10000e+01 % = 6 3 

With counter k set to 1 in step 250, a first derived system is found by performing 
step 252, thereby eliminating the variable x\ from equations 2 and 3. The coefficients 'e,y 
and the quantities x bi of the first derived system are as follows: 



JxBuUCeCl r Orrrl 


v ariuoies 


e[0][0] =+.10000e+01*a 


o 

e u 


e[0][l] =+.10000e+01 




e[0][2] =+.10000e+01 


°e\3 


b[0] =+.10000e+01*a+.20000e+01 


% 


e[l][0] =+.00000e+00 




e[l][l] = -.10000e+01+.10000e+01*a 


en 


e[ 1 ] [2] = -. 1 0000e+0 1 +. 1 0000e+0 1 *a 


e23 


b[l] =-.20000e+01+.20000e+01*a 


X b 2 


e[2][0] =+.00000e+00 


^31 


e[2] [ 1 ] = -. 1 0000e+0 1+. 1 0000e+0 1 *a 


l ^32 


e[2][2] = -.10000e+01-.10000e+01*a 


^33 


b[2] = -.20000e+01 


% 



The above first derived system for system Si 9 when written in normal algebraic 
form, appears as: 

ax\ + X2 + = a + 2 

(a -l)x 2 + (a - l)x 3 = 2(a-l) 

(a - l)x 2 -(a+ 1)jc 3 = -2 

By repeating steps 250 to 260, the method 200 calculates the second derived 
system for system Si as follows: 
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Reduced Form 



e[0 
e[0 
e[0 
b[0 
e[l 
e[l 
e[l 
b[l 
e[2 
e[2 
e[2 
b[2 



[0] =+.10000e+01*a 
[l] = +.10000e+01 
[2] = +.10000e+01 

= +. 1 0000e+0 1 *a+.20000e+0 1 
[0] = +.00000e+00 
[1] = -.10000e+01+.10000e+01*a 
[2] = -.10000e+01+.10000e+01*a 

= -.20000e+01+.20000e+01*a 
[0] = +.00000e+00 
[1] = +.00000e+00 

[2] = +.20000e+01*a-.20000e+01*a*a 
= +2.0000e+00*a-2.0000e+00*a*a 



Variables 

°en 
° e , 3 



e 2 \ 



x e 2 2 



e 2 3 
2 e 32 

2 ei 3 = l: 



33 



2 b 3 = r 3 



or alternatively 



ax\ + x 2 + X3 — a + 2 

(a -1)jc 2 + (a -l)x 3 = 2(a -1) 

-2a (a-l)x 3 = -2a(a- 1) 



Performing the back substitution step 280 the numerators n and the 
denominators /,-,■ can be found. In particular, from the last equation of the second derived 
system the numerator r 3 and the denominator 1 33 are as follows: 

1 33 = -2a (a - 1) and r 3 = -2a (a — 1). 

Substituting numerator r 3 and denominator / 33 into the second equation, we get: 
Reduced Form 

e[l][l] = -.20000e+01*a+.40000e+01*a*a-.20000e+01*a*a*a 
b[l] = -.20000e+01*a+.40000e+01*a*a-.20000e+01*a*a*a 

or 

l 22 = -2a (1 - 2a + a 2 ) and r 2 = -2a (1 - 2a + a 2 ). 



Variables 
hi 
r2 



l 
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In the final back substitution we get 

Reduced Form Variables 

e[0][0] = -.40000e+01*a*a*a+.12000e+02*a*a*a*a-.12000e+02*a*a*a*a*a 

+.40000e-K)l*a*a*a*a*a*a l n 
b[0] =-.40000e+01*a*a*a+.12000e+02*a*a*a*a-.12000e+02*a*a*a*a*a 

+.40000e+01*a*a*a*a*a*a r x 

producing thereby 

In = -4a 3 (1 - 3a + 3a 2 - a 3 ) and r x = -4a 3 (1 - 3a + 3a 2 - a 3 ). 



In a similar manner, the first derived system of system S2 may be written 

follows: 

Reduced Form Variables 



e[0] 


[0] =+.10000e+01*a 


°e u 


e[0] 


[1] =+.20000e+01 


%2 


e[o; 


[2] = +.00000e+00 


°e n 


b[0 


] = +.10000e+01*a+.20000e+01 


% 


e[i; 


[0] = +.00000e+00 


ei\ 


e[r 


[1] = -.40000e+01+.20000e+01*a 


en 


e[r 


[2] = +.00000e+00 


e23 


b[l 


] = -.40000e+01+.20000e+01*a 


X b 2 


e[2; 


[0] = +.00000e+00 


x e 3X 


e[2; 


|[1] = +.10000e+01*a 


^32 


e[2; 


|[2] = -.10000e+01*a 


l e 3 3 


b[2 


] = +.00000e+00 


l b 3 



ax\ + 2x2 = a + 2 
2 {a -2)x 2 = 2 {a -2) 
0x2 — 0x3 = 0 



The second derived system for system S2 is as follows: 
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Reduced Form 



Variables 



or 



e[o; 


][0] =+.10000e+01*a 


°e u 


e[o; 


][1] =+.20000e+01 


°e i2 


e[o; 


| [2] = +.00000e+00 


°en 


b[0 


] =+.10000e+01*a+.20000e+01 


% 


e[r 


|[0] =+.00000e+00 




e[i: 


|[1] = -.40000e+01+.20000e+01*a 


X e 22 


e[i; 


| [2] = +.00000e+00 










e[2; 


[0] = +.00000e+00 


2 en 


e[2j 


[1] = +.10000e+01*a 


2 

^32 


e[2] 


[2] = +.40000e+01*a-.20000e+01*a*a 


2 e 33 = / 33 


b[2! 


| =+.40000e+01*a-.20000e+01*a*a 


2 b 3 = r 3 


ax\ 


+2x2 — a + 2 




2(a 


- 2)x 2 = 2{a - 2) 





2a (2 - a)x 3 = 2a(2- a) 

Again performing the back substitution step 280 with system S2 the numerators n 
and the denominators /„ can be found. The numerator r 3 and the denominator / 33 are as 
follows: 

/ 3 3 = 2a (2 — a) and r 3 = 2a (2 — a). 

Substituting numerator r 3 and denominator / 33 into the second equation, we get: 



or 



Reduced Form 



JP9-1 999-0272 



Variables 



e[l][l] = -.16000e+02*a+.16000e+02*a*a-.40000e+01*a*a*a l 22 
b[l] = -.16000e+02*a+.16000e+02*a*a-.40000e+01*a*a*a r 2 

l 22 = -4a (4 - 4a + a 2 ) and r 2 = -4a (4 - 4a + a 2 ). 
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In the final back substitution we get 



Reduced Form 



Variables 



e[0][0] =-.64000e+02*a*a*a+.96000e+02*a*a*a*a-.48000e+02*a*a*a*a*a 



b[0] =-.64000e+02*a*a*a+.96000e+02*a*a*a*a-.48000e+02*a*a*a*a*a 

+.80000e+01*a*a*a*a*a*a r x 

or 

h 1 = -8a 3 (8 - 1 2a + 6a 2 - a 3 ) and n = -8a 3 (8 - 1 2a + 6a 2 - a 3 ). 

Performing step 290, the expressions * (r,) 2 and (/,,) 2 * (r,-)i are calculated 
and reduced to their reduced forms. For example, calculating (/ 2 2)i * (7-2)2 gives the 
following: 

(/22>i * (7-2)2= (-.20000e+01*a+.40000e+01*a*a-.20000e+01*a*a*a) * 



(-.16000e+02*a+.16000e+02*a*a-.40000e+01*a*a*a) 

= +.32000e+02*a*a-.32000e+02*a*a*a+.80000e+01*a*a*a*a 

-.64000e+02*a*a*a+.64000e+02*a*a*a*a-.16000e+02*a*a*a*a*a 
+.32000e+02*a*a*a*a-.32000e+02*a*a*a*a*a+.80000e+01*a*a*a*a*a*a 

= +.32000e+02*a*a-.96000e+02*a*a*a+.10400e+03*a*a*a*a 
-.48000e+02*a*a*a*a*a+.80000e+01*a*a*a*a*a*a 



= +.32000e+02*a*a-.64000e+02*a*a*a+.32000e+02*a*a*a*a 

-.32000e+02*a*a*a+.64000e+02*a*a*a*a-.32000e+02*a*a*a*a*a 
+.80000e+01*a*a*a*a-.16000e+02*a*a*a*a*a+.80000e+01*a*a*a*a*a*a 

= +.32000e+02*a*a-.96000e+02*a*a*a+.10400e+03*a*a*a*a 
-.48000e+02*a*a*a*a*a+.80000e+01*a*a*a*a*a*a 



+.80000e+01 *a*a*a*a*a* 



a 



Similarly, calculating (^2)2 * (^2)1 gives the following: 



(/22>2 * (^2)l 



(-.16000e+02*a+.16000e+02*a*a-.40000e+01*a*a*a) * 
(-.20000e+01 *a+.40000e+01 *a*a-.20000e+01 *a*a*a) 
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Step 290 similarly calculates the expressions (/,•/) i * (r,) 2 and (/„) 2 * (r,-)i for i=l 
and i=3. A simple string comparison of (/ 2 2)i * (r 2 ) 2 with (/ 22 ) 2 * (^2)1, performed in step 
300, shows that these expressions match. By repeating the comparison of (/,-,■) 1 * (V/) 2 with 
(I ah * (rj)\ for i=l and i=3 9 and finding that the expressions match for each i = 1,2 and 3, 
it can be shown that system £1 is equivalent to system £2. 

Embodiments of the invention can be implemented within compilers, for example. 
As is well known, a compiler generates machine executable object code from high-level 
source code, written in languages such as C++. 

The foregoing describes only some embodiments of the present invention, and 
modifications and/or changes can be made thereto without departing from the scope and 
spirit of the invention, the embodiments being illustrative and not restrictive. For 
example, the equivalence of more than two sets of simultaneous linear algebraic 
equations may be determined by pair-wise comparing the sets for equivalence. 
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